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Wc present the results of Monte-Carlo simulations of the charge carrier transport in a disordered molecular 
system containing spatial and energetic disorders using the dipolar glass model. Model parameters of the 
material were chosen to fit a typical polar organic photoconductor polycarbonate doped with 30% of aro- 
matic hydrazone, whose transport properties are well documented in literature. Simulated carrier mobility 
demonstrates a usual Poole-Frenkel field dependence and its slope is very close to the experimental value 
without using any adjustable parameter. At room temperature transients are universal with respect to the 
electric field and transport layer thickness. At the same time, carrier mobility does not depend on the layer 
thickness and transients develop a well-defined plateau where the current does not depend on time, thus 
demonstrating a non-dispersive transport regime. Tails of the transients decay as power law with the expo- 
nent close to -2. This particular feature indicates that transients are close to the boundary between dispersive 
and non-dispersive transport regimes. Shapes of the simulated transients are in very good agreement with 
the experimental ones. In summary, we provide a first verification of a self-consistency of the dipolar glass 
transport model, where major transport parameters, extracted from the experimental transport data, are 
then used in the transport simulation, and the resulting mobility field dependence and transients are in very 
good agreement with the initial experimental data. 



I. INTRODUCTION 

Dipolar glass (DC) model provides a natural and at- 
tractive description of major transport properties of polar 
amorphous organic materials. Initially, it was suggested 
for an explanation of the universal experimental observa- 
tion of the Poole-Frenkel (PF) mobility field dependence 



In /.J = const + 6S^/^ 



(1) 



in amorphous organic materials in a broad field range 
spanning almost two decades 10^ — 10^ V/cm.^^^ This 
model could be considered as a further development of 
the Gaussian Disorder Model (GDM),"* which assumes 
that charge transport in organic materials occures as 
a series of hops in random Gaussian energy landscape 
U (r ) , created by contributions from randomly packed 
molecules. Long ago Borsenberger and Bassler^ demon- 
strated that the GDM explains a strong polarity effect 
on carrier mobility by taking into account the dipolar 
as well as van der Waals energetic disorder. Yet, by the 
very nature of the GDM, it does not incorporate the most 
striking feature of the dipolar disorder, i.e. a long range 
spatial correlation of the random energy landscape,^''' 
and for this reason cannot adequately explains emergence 
of the PF dependence. In contrast with the GDM, DG 
model explicitly takes into account spatial correlations 
and naturally provides the PF field dependence in a wide 
field range. """^^ 



Initially, the DG model has been dubbed "Correlated 
Disorder Model" (CDM)^, stressing the most important 
qualitative difference of the model with the GDM, but 
the proposed name describes better the true nature of 
the model. Indeed, as it was shown in Ref. 7, different 
spatial correlations of the random energy landscape lead 
to the different field dependences of the mobility, hence 
the CDM name is too vague. The particular model, sug- 
gested and studied in Refs. 7 and 9, is best suited for 
the description of transport properties of polar organic 
materials containing molecules having permanent dipole 
moments. 

For the study of hopping charge transport in organic 
glasses a lattice version of the DG model is usually con- 
sidered with sites of the regular lattice occupied by ran- 
domly oriented dipoles, and some of sites are considered 
as transport ones. Earlier transport simulations for the 
DG model have been done for the case of the totally filled 
lattice where the fraction c of sites occupied by transport 
molecules was equal to 1. Such simulations suggested 
that the temperature and field dependence of the mobil- 
ity may be described by the phenomcnological equation 
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with Co ~ 0.78 and F w 2.^ Here cr is the magnitude of 
the energetic dipolar disorder and a is the lattice scale. 
This particular relation has been widely used for a de- 
scription of transport properties of various materials, we 
mention here just some recent references. ^^^^^ In some 
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papers a comparative analysis of the data using both the 
GDM and DG model was carried out and it was found 
that the DG model provides more coherent description of 
charge transport. Still, a proper experimental test 
of the DG model is far from completed. 

First, at the moment no simulation has been carried 
out for c < 1. The lattice model with c = 1 is ad- 
equate for description of low molecular weight organic 
glasses, while polymeric materials, especially molecularly 
doped polymers, are certainly better described by the 
case c < 1. In this paper we provide results of the trans- 
port simulation for a partially filled lattice with c = 0.3 
to verify the validity of the PF relation for that case. 
Numerous experiments demonstrate that the PF rela- 
tion is generally valid for c ~ 0.2 — 1; in addition, there 
is no theoretical reason to doubt its validity. Indeed, a 
general picture of the carrier hopping in the correlated 
energy landscape suggests that major features of quasi- 
equilibrium non-dispersive transport is governed by the 
critical valleys with the size^° 
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If the typical distance between transport sites ~ a/c^^^ 
is much less than 
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one should not expect any significant modification of the 
PF dependence. It is worth noting that according to 
this line of reasoning, with the decrease of c the mobility 
dependence initially begins to deviate from the strict PF 
law in the high field region. For w 1 x 10^ V/cm, 
a ~ 0.1 eV, and a ~ 1 nm inequality (4) gives c ^ 0.1 
at room temperature. Nonetheless, a direct check of the 
validity of relation ( 1 ) for a partially filled lattice is highly 
desirable. 

Next, the most important reason to perform the sim- 
ulation for c < 1 that uses set of transport parameters 
extracted from experimental data for a particular trans- 
port material is to provide a test of the self-consistency 
of the transport model. For quite a long time the GDM 
and DG model are used for the description of transport 
properties of organic materials, i.e. relevant transport 
parameters (mostly a) are calculated using experimen- 
tal mobility data according to the prescriptions of the 
particular model. All such prescriptions are provided by 
the phenomenological relations, based on the extensive 
Monte Carlo simulations."*'^ In this situation a very nat- 
ural suggestion should be to perform a simple but very 
convincing test for the self-consistency of the model: to 
carry out the simulation for a given set of parameters (ex- 
perimental ones) and then compare the results (mostly, 
the simulated mobility field dependence) with the ex- 
perimental data. To the best of our knowledge, such 
comparison has been done neither for the GDM nor DG 
model. In this paper we present the result of such test 



for the particular set of transport parameters, extracted 
from the experimental mobility data for a typical organic 
semiconducting glass containing 30% (by weight) hole 
transporting aromatic hydrazone DEH in polycarbonate 
(PC) matrix, and then compared our results with the ex- 
perimental ones.^^^^'^ Experimental results, provided in 
different papers, agree well, thus justifying our choice of 
the reference transport material. 

Last, earlier Monte-Carlo simulations^^ demonstrated 
that time-of-flight (TOF) current transients, predicted 
by the DG model, reproduce very well the experimental 
transients: there is a short spike reflecting initial spa- 
tial and energetic relaxation of carriers, followed by a 
flat plateau signaling emergence of the quasi-equilibrium 
transport with the mean carrier velocity independent of 
time, and then a final tail marking the arrival of carriers 
to the collecting electrode and typically having algebraic 
form I{t) oc l/t^ with /3 ~ 2.0 - 2.5. 

At the same time, a general behavior of transients was 
not well studied because the paper^^ was entirely devoted 
to the study of an important but very particular effect of 
deep traps on the shape of transients. In this paper we 
study TOF current shapes predicted by the DG model 
via Monte Carlo simulation for a partially filled lattice 
and compare simulated transients with the experimental 
ones for 30% DEH:PC. Again, we use this comparison as 
a test of the self-consistency of the DG model. 

Realization of this program could provide the DG 
model with more strong footing. 



II. BASICS OF THE DG MODEL AND SIMULATION 
DETAILS 

Major assumptions of the DG model are well 
known. ^'^'*°'^* Usually the model is considered as a reg- 
ular cubic lattice with sites occupied by randomly ori- 
ented dipoles having dipole moment p. We assume that 
orientation of dipoles is static and not affected by the 
applied uniform electric field E. The dipoles generate 
a random spatial distribution of the electrostatic po- 
tential Lp{r) which immediately translates to the ran- 
dom site energy U{r) = eip{r) with zero mean and rms 

1/2 

a = (U'^{r)) = 0.05-0.1 eV, depending on a andp.^-'^ 
Particular energy distribution was obtained by the sum- 
mation of contributions of all dipoles using the Ewald 
method^^ and assuming the periodic boundary condi- 
tions for the distribution of dipoles in the basic sample; 
we assumed also that every site is occupied by a dipole. 
If the fraction of sites, occupied by dipoles, is not very 
low and the average distance between neighbor dipoles 
is comparable to a, then the resulting distribution of U 
has a Gaussian form; for low concentration of dipoles 
the tail of distribution decays more slowly.^'^^'^*" If the 
energy landscape has a Gaussian distribution, then the 
only relevant parameter is cr; fraction of sites, occupied 
by dipoles, the number of different sorts of dipoles, and 
the very nature of dipoles (i.e., are they transport sites 



Charge transport in doped polycarbonate 



3 



or not) are irrelevant for the problem. 

Random energy landscape provided by dipoles has a 
particular property which is in drastic contrast to the 
GDM: a binary correlation function C(r) = {U{'r)U{0)) 
decays slowly with distance, C(r) cx a/r, r 3> a,^'^ 
while for the GDM this function is zero everywhere apart 
from r = 0, where C(0) = cr^. Long range corre- 
lation naturally leads to the PF dependence: the ID 
transport model predicts that if C(r) oc l/r", then 
In^i a £:"/("+!). 7 

In typical organic polar materials the actual random 
energy landscape U{f) incorporates contributions from 
other sources of energetic disorder: quadrupolar,^*'^^ 
van-der-Waals,'^° etc. The dipolar contribution stands 
apart because its correlation function decays more 
slowly than corresponding functions for other sources of 
disorder. ^^'■^^ This means that the dipolar disorder, when 
not negligible, provides the dominant contribution to the 
mobility field dependence for low and moderate fields. ^'■^^ 
This fact justifies the use of a pure dipolar disorder (i.e., 
the DC model) for description of the mobility field de- 
pendence in polar amorphous organic materials for a field 
range 10^ - 10^ V/cm. 

Our simulation technique is very close to the one, de- 
scribed previously.^ The only difference is that the frac- 
tion of sites, occupied by transport molecules is not c = 1 
but c = 0.3, thus providing connection with the ex- 
perimental data for a particular well-characterized sys- 
tem, 30% DEH:PC. ^^"22,33 ^gg^ ^^le simplest assump- 
tion that sites of the lattice are occupied by transport 
molecules independently with probability c; the use of 
the simplest hypothesis is mostly related to the lack of 
relevant experimental information on the local structure 
of organic glasses. 

In our simulations we examine a temporal dependence 
of the carrier velocity v(t), averaged over many carriers 
for a particular realization of the random energy land- 
scape U{r), and then over realizations of U{r). Current 
transient I{t) in the TOF experiment is directly propor- 
tional to v{t). Another important transport parameter, 
used for the calculation of the Monte Carlo carrier mo- 
bility, is a mean carrier velocity {v) = (i/tdrift), where 
idrift means the time for a carrier to reach the collecting 
electrode. 

We used the Miller- Abrahams (MA) hopping rate,'^''^ 
where the rate of transition from site i to site j is given 

by 
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where vq is the prefactor frequency, 



Uj - C/j < 
(5) 
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7 is a wave function decay parameter for transport sites; 
in organic materials 7a ~ 5 — 10.^'"^^ In the presence of 
applied electric field E the site energy Ui includes an ad- 
ditional term —efiE. It was found previously that details 
of the hopping rate are irrelevant to the emergence of the 
PF dependence (1), at least for the case c = 1.^ 



Relevant physical parameters for the simulation were 
taken from the experimental transport data for 30 wt% 
DEH:PC glass at room temperature. ^^"^^ We used g = 
0.13 eV, a = 7.7A, 27a = 11.8, and kT = 0.0252 eV. In 
these papers cr was calculated using the GDM analysis of 
the low field mobility. The corresponding DC value ctdg 
could be obtained multiplying by the correction factor 
10/9."*'^ Wc did not perform this correction because it is 
well known that the coefficient A in the relation 



ln(^/pto) \e^o = -A{a/kT) 



(6) 



valid for both GDM and DG model, is not a constant but 
depends on ^a?^ This dependence is rather weak around 
7a 5 — 10, but, nonetheless, it provides some uncer- 
tainty for a value of a, calculated from the experimental 
mobility data. In addition, there may be a weak depen- 
dence of A on c as well. In this situation it is difficult 
to make a reliable but rather small correction from the 
GDM to DG model. 

Exceptions are Figs. 1 and 3, where simulation was 
performed using cTcorr = 10/9 cr = 0.144 eV as well, in 
order to estimate the agreement between our simulation 
and experimental mobility data for 30% DEH:PC. Sim- 
ulation time becomes prohibitively long to obtain good 
quality transients (especially long-time tails) using (Tcorr, 
while the corresponding simulation of a more robust inte- 
gral transport parameter (i.e., the mobility) still remains 
feasible. Nonetheless, there is no doubt that the major 
features of transients (such as universality, see Section 
IIIB) remain intact for the case of corrected cr. 

Quite frequently a polaron (or Marcus) hopping rate 
is considered as a more realistic alternative to the MA 
rate 



37,38 



exp 



{Uj -U, + 2Etf 
SEbkT 



(7) 



where Et is the polaron binding energy and we omit 
here the dependence on distance. For typical organic 
molecules used in transport materials intramolecular con- 
tribution to Eb is 100-150 mcV,"^^ and comparable inter- 
molecular contribution could be expected. In this sit- 
uation a major difference with the MA rate is an addi- 
tional contribution £'f,/2 to the effective activation energy 
2Aa^ /kT. Significant effect on the shape of transients 
(as well as on the mobility field dependence) could be 
expected for rather exotic case of small Et ~ 10 — 30 
meV, where transport is dominated by the quadratic 
term {Uj — Ui)'^ in Eq. (7) and occurs in the inverted 
regime. ■^^ This statement is supported by our prelimi- 
nary data. Detailed study of transients for the polaron 
hopping rate will be published elsewhere. At the same 
time, a possible polaronic contribution to the transport 
activation energy provides an uncertainty for the value 
of cr, estimated from the experimental data. Unfortu- 
nately, at the moment there is no reliable method to 
extract polaronic contribution to the mobility temper- 
ature dependence. For this reason we assumed that all 
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activation dependence originates from the random energy 
landscape. 

In order to check the sensitivity of the simulation data 
to a finite size effect, mentioned in Refs. 40 and 41, 
we performed simulation for basic samples with size 
equal to 50^, 100"^, and 150'^ lattice sites and found that 
the parameter, most sensitive to A^, is the absolute value 
of the mobility. Mobility field dependence and relative 
shape of the transients are much less sensitive. For ex- 
ample, slope of the mobility field dependence S = 2^tj2 
in Fig. 1 varies by 5% only and general features of tran- 
sients remain the same. In addition, the finite size effect 
has been demonstrated for i? = as a deviation of the 
mean carrier energy en from the value £oo = — cr^/fcT, 
which is a mean carrier energy for E = and infinite 
number of sites in the case of the Gaussian distribution; 
for a finite N en > £oo- The reason for the deviation 
is scarcity of sites with low energy. However, for E > 
mean carrier energy increases with E and becomes more 
close to the maximum of the density of states. For this 
reason the finite size effect becomes less severe for the 
simulation with E > 0. 

We should note two earlier papers where the Monte 
Carlo simulation of the hopping transport in amorphous 
organic materials for c < 1 has been carried out assum- 
ing random distribution of hopping sites. ^'^'^^ Contrary 
to our approach, simulation in Ref. 43 has been carried 
out for the GDM. Sin and Soos'*"' did not simulate tran- 
sients and considered only the case of moderate disorder 
cr/kT < 3, while for 30% DEII:PC and room temperature 
a/kT > 5. In addition, they generated the correlated en- 
ergy landscape in somewhat artificial way according to 
the method suggested by Gartstcin and ConwcU."'^ This 
method gives short range spatial correlation and cannot 
reproduce the proper spatial correlation in polar organic 
materials. 

Some earlier papers considered the concentration de- 
pendence of the hopping mobility in organic materials us- 
ing various analytic approaches. Unfortunately, such 
approaches do not provide any information on the shape 
of transients. 
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FIG. 1. Comparison of the experimental mobility field de- 
pendence (denoted by □, data are borrowed from Refs. 
21, 23, and 33; we do not discriminate points from differ- 
ent papers, they are agree well with each other) and simu- 
lated ones (DG model for a — 0.13 eV (T) and corrected 
aeon = 0.144 eV (■), and for the GDM (•), correspondingly). 
The topmost curve (a) shows the simulation data for the 
equivalent 100% filled lattice with lattice parameter equals 
to a/0.3^^''. Plotted is the logarithm of the ratio n/ns, where 
fj,5 is the corresponding mobility for E — 5 Y/fim (matching 
point is indicated by the arrow). Hence, the absolute value of 
the mobility is irrelevant, and only the slope of the mobility 
field dependence is compared. Note that in this figure there 
are no fitting parameters. Dotted lines are used as guides for 
an eye. Inset shows our mobility curve for acorr = 0.144 eV 
(hence, a/kT = 5.7) and the mobility curve for a/kT = 5.8 
(o), borrowed from Fig. 3 of Ref. 50; for i? = 5 V/ixm in our 
case (eaE/a)^^^ = 0.163 and for the reference mobility value 
for data from Ref. 50 we used, instead of fi^, the mobility, 
simulated for the same value of (eaE /a)^^^ . 



III. RESULTS AND DISCUSSION 

A. Mobility field dependence 

We found that the PF field dependence is valid for c = 
0.3 as well as for a fully filled lattice (Fig. 1). Deviation 
from the PF dependence begins around 1 x 10^ V/cm, re- 
fiecting the specific feature of the Miller- Abrahams hop- 
ping rate^ and effect of the decrease of the concentration 
of transport sites, mentioned in the Introduction. This 
figure also shows that the DG model successfully passes 
the self-consistency test and simulation with parameters, 
extracted from experimental data, very well reproduce 
initial experimental field dependence. Simulation using 
(Tcorr docs not differ significantly from the simulation us- 



We also carried out simulation for the equivalent to- 
tally filled lattice, i.e. the lattice with Ce = 1 and 
ttg = a/c^l'^ K, 1.5a which simulates the same material 
30% DEH:PC in the traditional lattice gas model. ^ If we 
compare average velocities, then for the partially filled 
lattice we have v = va/r, while for the equivalent lattice 
(for the same physical thickness of the transport layer 
and electric field E) — v^ae/Te, where raw output of 
the Monte Carlo simulation v is the dimensionless veloc- 
ity, measured in the units of a/r (here t = i'q^ exp (27a) 
is the characteristic hopping timcscalc). For the equiv- 
alent lattice Te = exp {2jae), and the ratio of merit 
is 
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If the approximation of the equivalent lattice is valid, 
then i? ~ 1. In our case, for E = 10 V/fim we obtain 
R = 1.85; this is a reasonably good agreement (note, 
that Ve/v « 230). However, for higher fields the agree- 
ment becomes consistently poorer (see the topmost curve 
in Fig. 1). At the same time, almost all discussions of 
the comparison of transport properties of organic mate- 
rials with different concentrations of transport sites ex- 
tensively use, albeit sometimes implicitly, the equivalent 
lattice approximation (see, for example, recent Refs. 32 
and 49). Conclusions of such discussions should be reex- 
amined. 

It is natural to assume that the mobility curve for the 
totally filled lattice could be made more close to the ex- 
perimental one by the addition of the effective spatial 
disorder according to the prescription of Ref. 4 because 
the addition of such disorder leads to the decrease of the 
slope S of the mobility field dependence. However, in 
this approach the spatial disorder is introduced in a pure 
phenomenological manner without any possibility to pre- 
dict a magnitude of disorder, necessary for the reproduc- 
tion of the proper mobility field dependence. Quite the 
contrary, in our approach a simple use of the proper con- 
centration of transport molecules immediately provides a 
good description of the mobility field dependence. 

Our results also provide good illustration on the impor- 
tance of long range spatial correlations for the develop- 
ment of the proper mobility field dependence. Sometimes 
in the literature one can find a statement that short range 
correlations are quite sufficient for the development of the 
field dependence which is almost indistinguishable from 
the true PF dependence in a wide field range and, hence, 
for a description of transport properties of organic mate- 
rials (a recent example is provided in Ref. 50) In order to 
understand the importance of long range correlations it 
is useful to compare our Fig. 1 (inset) and Fig. 3 of the 
cited paper. In our case the PF dependence starts imme- 
diately from the lowest tested field, while for the short 
range disorder the slope for the eaE/a < 0.3 is almost 
zero, or even negative, though, probably, the reason for 
this unusual behavior is small thickness of the simulated 
transport layer. For correlations that are negligible for 
distances greater than 3a (Fig. 2 of Ref. 50) the result- 
ing mobility field dependence has to be very close to the 
GDM dependence, and this is clearly seen in Fig. 1. 
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FIG. 2. Universality of I{t) (or temporal dependence of v): 
for electric field in the range 5 - 200 Y/fim and L — 20 000 a, 
here different lines correspond to the different values of E (a); 
and for layer thickness in the range 2 000 a - 20 000 a, E = 5 
Y/fim, here different lines correspond to the different values 
of L (b). For both plots kT/a = 0.19 (room temperature). 



B. Current transients 

It was already established previously that in our case 
charge transport occurs in the non-dispcrsivc regime; for 
this reason mobilities calculated in different ways (as a 
true Monte Carlo mobility ^ = (v) /E calculated by mean 
carrier velocity (v), or using two experimental methods 
of the mobility calculation, estimating either the time of 
intersection of asymptotes to the plateau and tail of the 
transient, or the time of the transient to reach one half of 
the plateau value), are very close. Taking into account 



all relevant parameters for a non-dispersive regime with 
idrift oc L (here L is a thickness of the transport layer) , a 
general shape of the transient may be written as 

vit) = (v) F3 i{v) t/L, kT/a, E/Es) , (9) 

i.e., as a three-parameters scaling function, where Es is 
some characteristic field. We found that at room tem- 
perature {kT/a ~ 0.19) transients demonstrate a univer- 
sality with respect to L and E (see Figs. 2-4), found 
previously in experiments. ^^'^^'^^ In agreement with ex- 
periment, transients are not universal with respect to T. 
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FIG. 3. Universality of v{t) for electric field in the range 
55 - 200 Y/^im and L = 20 000 a, here different lines corre- 
spond to the different values of E; transient are simulated for 
o"corr = 0.144 eV. Inset shows transients for E = 55 V//im, 
simulated for the GDM (broken line, a = 0.13 eV) and DG 
model (solid line, aeon ~ 0.144 eV); ts and Vs provide a suit- 
able normalization. 

Universality was also confirmed in a more limited field 
range for the transients simulated using CTcoit = 0.144 
eV (Fig. 3). In this case transients are more dispersive, 
yet still less dispersive in comparison to the correspond- 
ing GDM transients, calculated for a = 0.13 eV. Fig. 4 
shows the field dependence of the parameter 

^^h/l^^ (10) 

tl/2 

which provides a simplest robust integral characteristic of 
the transient shape. Here to is the time of the intersection 
of the asymptotes to the plateau and tail of the transient, 
and ti/2 is the time for the current to decay to one half 
of the plateau value. If Eg does not depend on L (i.e., 
if it is some microscopic field), then Eq. (9) implies the 
universality with respect to L, in agreement with the 
simulation result, but the universality with respect to E 
means in addition that the third parameter in Eq. (9) is 
irrelevant, and the two-parametric scaling takes place 

v{t)^{v)F2{{v)t/L,kT/a), (11) 

where the dependence of v{t) on E is exclusively provided 
through {v). 

It is not clear, why the third parameter in Eq. (9) is ir- 
relevant. Considering the nature of the problem, we may 
suggest two possible candidates for Eg: E^^^ = a /ere 
and = e/r^. In fact, if we use Eg'^^ in Eq. (9), 

(2) 

then E/Es docs not depend on E, and this is exactly 



FIG. 4. Dependence of the parameter W on E for room tem- 
perature; inset shows the dependence of VK on T for i? = 10 
V/nm. 

the universal behavior of Eq. (11), but, again, why one 
should use i?s^'* instead of iJs^' is unclear. 

Typically, the universality of transients is considered 
as an inherent property of the dispersive transport. In 
our case validity of Eq. (11) means that the mean carrier 
velocity (and, hence, the mobility) does not depend on L 
and, in this respect, the transport is also non-dispersive. 

The most interesting feature is the behavior of the tail 
of the transients. We found that at room temperature 
the tail obeys the law 

v{t)(xl/t^, (12) 

with parameter /3 being very close to 2 (Fig. 5). It is 
quite probable that (3 is even more close to 2 because de- 
termination of the tail parameters is very sensitive to in- 
evitable statistical errors. For the equivalent 100% filled 
lattice /3 strongly depends on E and transients arc not 
universal (Fig. 5b, inset). 

Close values of (3 are observed for other organic glasses, 
too (see inset in Fig. 5a for the transient in fiuorene- 
arylamine copolymer at room temperature and E = 
1.9 X 10^ V/cm).^"* For this particular glass the univer- 
sality with respect to E takes place, too, and the tran- 
sients demonstrate well-defined plateau indicating the 
non-dispersive transport. In Ref. 52 it was found that 
/3 « 3 and again is independent of E and L. Note, 
though, that the the actual value of (3 is sensitive to the 
time range, available for the analysis. Typically, in ex- 
periments the tested time range is not very wide and the 
transients decay just by one order of magnitude (see in- 
set in Fig. 5a) . In this situation true values of /3 for very 
long times might be even more close to 2. 
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FIG. 5. Typical shape of the transient for _E = 10 V//xm, L = 
20 000 a, inset shows the transient for the fluorene-arylamine 
copolymer, experimental data is borrowed from Fig. 2 in Ref. 
54 (a), and the dependence of the parameter /3 in Eq. (12) on 
the electric field, inset shows the corresponding dependence 
for 100% filled equivalent lattice (b). 



It is worth noting that the power law decay of tran- 
sients has been observed even in liquid crystals, and again 
it was preceded by a well-defined plateau, though in 
that case /? « 3 or 4, depending on the particular liquid 
crystalline phase. Hence, even in much more locally or- 
dered materials the tail of the transient does not follow 
the exponential law, inherent for classical diffusion. 

Power law decay of the transient is typically associated 



FIG. 6. Simulated transient for iJ = 36 V//im and L = 
20 000 a (solid line) and transients 3 and 4 from Fig. 1 in Ref. 
23 (broken and dotted lines, correspondingly); both experi- 
mental transients have been measured at i5 = 40 V/^m but 
for slightly different thickness of the transport layer (14/im 
and 15/im) . Experimental transients were re-scaled to obtain 
the best fit. Inset shows the more dispersive GDM transient 
with no visible plateau for the same E. 



with dispersive transport, where 



I{t) cx 



t < tT, 



t-(l+"), t>tT, 



(13) 



where is some characteristic transport time and < 
a < 1. If /3 « 2, then a = /3 - 1 w 1. Formally, this 
means that wc have the boundary case between dispersive 
and non-dispcrsivc transport, and the carrier mobility 
docs not depend on the layer thickness L, because fi oc 
j^i-i/a — jji — const. Hence, the power law decay of the 
transients with (5 k, 2 and independence of {v) on L are 
closely related and mutually consistent phenomena. 

We checked an agreement between shapes of the simu- 
lated transients and experimental ones, measured in 30% 
DEH:PC layers. ^■^ Result of the comparison is shown in 
Fig. 6. Difference between values of E (and, hence, an 
average carrier velocity) and L in experiment and simu- 
lation {L = 15.4 ^m) is not important due to transient 
universality: in this procedure wc compare only the rel- 
ative shapes of the transients. 

The GDM shows poorer universality in the same field 
and thickness range (Fig. 7). Actually, for most values 
of E the GDM transients are so dispersive that reliable 
determination of the transit time in double linear coordi- 
nates and, hence, calculation of W is not possible. This is 
a very typical situation: if wc extract the relevant trans- 
port parameters for the GDM from the experimental data 



Charge transport in doped polycarbonate 



8 





0.5 1 1.5 2 2.5 3 

<v>t/L 



FIG. 8. Dependence of the mean carrier velocity (v) on the 
thickness L (number of lattice planes) of the transport layer 
for E ^ 55 V/fim in the GDM (o) and DG model (□), cor- 
respondingly. Velocity is normalized by the corresponding 
velocity (w)2oooo 20 000 lattice planes. Lines are provided 
as guides for an eye. 



cause in both cases the corresponding values of a are not 
equal for t < tx and t > tx- This fact was already men- 
tioned by Bassler.^ Fig. 8 provides additional support 
to unusual dispersion of transients in organic glasses, be- 
cause if Eq. (13) is valid, then {v{L)) oc L^~^/°', which 
docs not agree with Fig. 8. 



IV. CONCLUSION 



FIG. 7. Check of the universality for the GDM transients: 
transients for electric field in the range 5 — 200 V / and 
L = 20 000 a (a) ; transients for layer thickness in the range 
2 000 a - 20 000 a and E = 36 V/fim (b). 



and then use these parameters for the simulation of the 
GDM transients, the output of the simulation produces 
much more dispersive transients than the experimental 
ones.^^ Dispersive regime for the GDM transport can be 
illustrated by the dependence of the mean carrier veloc- 
ity (v) on L (sec Fig. 8). Mean velocity monotonously 
decreases with L for the GDM, while for the DG model 
it is almost a constant for L > 2 000 a. This observation 
supports the conception of the breakdown of the this par- 
ticular type of universality for the GDM, because for the 
validity of Eq. (11) the mean velocity should be indepen- 
dent of L. At the same time, dispersive transport in the 
GDM and DG model is not described by Eq. (13), be- 



We carried out Monte Garlo simulation of the hopping 
charge transport in the DG model with parameters (a 
and 7a) directly taken from the experimental data for 
the archetypical molecular doped polymer 30% DEH:PC 
and assuming the simplest independent random distribu- 
tion of DEH transport molecules at the sites of a cubic 
lattice. We found that the mobility field dependence for 
moderate field in the DG model for c = 0.3 retains its 
PF form, in good agreement with the general picture of 
the hopping transport in the correlated energy landscape. 
The most important parameter, namely slope of the mo- 
bility field dependence, was found to be very close to the 
experimental one. 

Hence, for the DG model the simulation that uses the 
experimental parameters, derived from the temperature 
dependence of the low-field mobility (a) and dependence 
of low field mobility on concentration of transport sites 
(7a), is capable to reproduce well the mobility field de- 
pendence and shape of the transients. This agreement 
becomes even more convincing if one takes into account 
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that for the model of equivalent 100% filled lattice the 
slope of the mobility field dependence and shape of tran- 
sient are in striking disagreement with the experiment 
(see Figs. 1 and 5). And then, the simplest natural mod- 
ification of the model (i.e., the use of partially filled lat- 
tice and assumption of the random independent spatial 
distribution of transport sites) immediately brings simu- 
lated curves very close to the experimental ones. There 
is no need to introduce additional spatial disorder apart 
from the random distribution of transport site. 

We would like to emphasize an importance of this fact 
because recent survey suggested that the origin of the 
energetic disorder in molecular doped polymers and very 
nature of the activation energy of the hopping transport 
in these materials should be reconsidered.^^ According to 
this study, activation energy of the hopping transport in 
molecularly doped polymers is mostly of the intramolecu- 
lar origin. At the same time, there is a general belief that 
the energetic disorder with the magnitude cr ~ 0.1 eV is 
typical for amorphous organic materials. Moreover, long 
range spatial correlation of the energy landscape is, at the 
moment, the only possible candidate for the explanation 
of the PF mobility field dependence. ""'^^ Our simulation 
shows that it is possible to reproduce a typical field de- 
pendence in polar organic material (i.e., the dipolar glass 
with a « 0.13 eV) without using adjustable parameters 
and taking only a simplest assumption of the random dis- 
tribution of hopping sites in the bulk of organic material. 
In this simulation the magnitude of the slope of the mo- 
bility dependence is directly dictated by the magnitude of 
the dipolar disorder. Hence, the good agreement between 
our simulation and the experiment^^'^'^''^'^ provides an ad- 
dition support in favor of the existence of correlated en- 
ergetic disorder with a ~ 0.1 eV in organic materials. At 
the same time, as it was stressed in Ref. 49, in polymers, 
doped with highly polar transport dopants, the magni- 
tude of the total disorder, which presumably includes a 
significant and variable dipolar contribution, and calcu- 
lated using the temperature dependence of the low field 
mobility, seemingly does not depend on the dopant con- 
centration. This disagreement, why the dipolar disorder 
clearly reveals itself in the mobility field dependence and 
yet its contribution to the mobility temperature depen- 
dence is hidden, poses a major puzzle for the problem of 
carrier transport in amorphous organic materials. More- 
over, recent study casts doubt on the validity of the use of 
correlated models for description of the charge transport 
in polymers, in contrast with the transport properties of 
low molecular weight materials. The analysis was based 
on rather indirect comparison of the measured current- 
voltage curves with the predicted ones. At the same 
time, much more simple and direct TOF experiments un- 
ambiguously demonstrate that mobility field dependence 
agrees with predictions of the correlated model in both 
polymers and low molecular weight glasses, without any 
qualitative difference.^ Quite possibly, that the disagree- 
ment found in Ref. 58, also stems from the implicit use 
of the lattice gas model. New experiments on the concen- 



tration dependence of transport parameters with careful 
analysis of the transients along the lines suggested in Ref. 
32, as well as more thorough study of the simulated mo- 
bility dependence on the concentration of transport sites, 
are extremely desirable. At the same time, further com- 
parison of the simulation data with experimental ones 
could be considered as a logical extension of our study. 
Possible tests may include attempts to describe trans- 
formation of transients and variation of the slope of the 
mobility field dependence with temperature and concen- 
tration of transport molecules. A major difficulty for the 
realization of that programme is the scarcity of raw tran- 
sient data in the literature. '^^ 

We showed also that at the room temperature tran- 
sients demonstrate universality with respect to L and 
and are non-dispersive according to the usual definition 
(carrier mobility does not depend on i, and the transient 
demonstrates a flat plateau if plotted in double linear / 
vs t coordinates). Yet the shape of the transients dif- 
fers significantly from the result of the simple diffusive 
approximation, which was suggested for the description 
of experimental TOF transients. ^^"^^ Two major differ- 
ences are the universality of the transients with respect 
to L and power law decay of the tail. It is worth not- 
ing that universality with respect to L naturally arises 
in the model where charge transport is described by the 
broad distribution of the effective carrier velocities, ^^'^^ 
but universality with respect to E cannot be explained by 
this approach. It is also worth to note that the develop- 
ment of a flat plateau and simultaneous power law tail of 
the current with /3 ~ 2 cannot be described by the mul- 
tiple trapping model, frequently used for the description 
of the hopping transport in random media. ^'^ 
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